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The two-dimensional J -J' dimerized quantum Heisenberg model is studied on the square lattice by 
means of (stochastic series expansion) quantum Monte Carlo simulations as a function of the coupling ratio 
a — J' /J. The critical point of the order-disorder quantum phase transition in the J -J' model is determined as 
Qc = 2.5196(2) by finite-size scaling for up to approximately 10 000 quantum spins. By comparing six dimer- 
ized models we show, contrary to the current belief, that the critical exponents of the J -J' model are not in 
agreement with the three-dimensional classical Heisenberg universality class. This lends support to the notion 
of nontrivial critical excitations at the quantum critical point. 
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Dimerized quantum spin systems are important examples of 
low-dimensional antiferromagnets featuring a quantum phase 
transition (QPT) fill which destroys a Neel ordered state by 
competition between different interactions. In contrast to 
other examples showing such criticality, in this class of mod- 
els the actual transition is triggered by nonisotropic couplings 
where the dimers 12811 are explicitly placed on the lattice. 
Because of the discovery of Bose-Einstein condensation of 
magnons in a magnetic field much effort has been spent to 
study their physics |01. 

The characteristics of the QPT in two-dimensional (2D) 
dimerized models have been investigated in detail. By map- 
ping to a nonlinear sigma model (NLSM) Istl it was argued 
that the transition is well described by the Heisenberg 0(3) 
classical universality class in three dimensions (3D). The role 
of Berry phase terms, which are present in the mapping to the 
NLSM, is argued to be irrelevant 101 and there are numerous 
numerical studies which support this claim. Examples include 
the CaVO lattice |5i], bilayer models and the 2D coupled 
ladder system 101 ■ 

Recently, the idea of deconfined quantum critical points 
has been put forward by Senthil et al. iQl who argue that 
there are, however, also important examples of QPTs where 
Berry phases and nontrivial excitations at the quantum criti- 
cal point can change the critical behavior These arguments 
are based on Heisenberg models with isotropic interactions 
exhibiting a transition between two ordered states such as an 
antiferromagnetic and a valence-bond solid phase, and numer- 
ical evidence for such a case was recently claimed by Sandvik 
J^. This idea that challenges the standard Landau-Ginzburg- 
Wilson framework of phase transitions was also found to be 
relevant in other systems such as classical dimers ifioll in 3D 
and has prompted further theoretical and numerical efforts 
some of which extend to different sce- 
narios or show that the field is still highly controversial. 

In this context, it is in any case somewhat surprising that 
also a specific 2D dimerized spin model, which we refer to as 



the J-,J' model, with nonisotropic interactions was suggested 
as a candidate for deconfinement at the quantum critical point 
by Yoshuika et al. flill (see also Ref. IUtIi ). In consequence 
this idea could lead to critical exponents characterizing the 
phase transition, which differ from those of the Heisenberg 
universality class in 3D. This conclusion was, however, ques- 
tioned llSll because of the close relation of the J -J' model to 
the ladder model. In order to resolve this conflict and to give 
arguments in favour of one or the other alternative we report 
in this Letter on quantum Monte Carlo (QMC) simulations of 
various dimerized models, which signal the emergence of an 
unconventional phase transition for the J- J' model. 

The ,J-J' model is defined on a square lattice with N = L"^ 
spins by the Hamiltonian 



TC — J ^ ^ SiSj + J' ^ ^ SiSj 



(1) 



Here, = (1/2) (o-^, ay, a^) denotes the usual spin-1/2 op- 
erator at lattice site i, and J and J' are the antiferromagnetic 
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Figure 1: (color online), (a) Visualization of the J- J' model on the 
2D square lattice. The quantum spin (5 = 1 /2) degrees of freedom 
live on a square lattice with different nearest neighbor couplings J 
and J' (thin and thick), (b) Similar for the plaquette model, favoring 
quadrumer formation. 
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Figure 2: (color online), (a) The Binder parameter Q2 and (b) the correlation length (,y/L for various lattice sizes from L = 8 to L — 72. (c) 
Scaling of the crossing points can be used to extract the critical coupling Oc. All quantities seem to converge to the same estimate. 



coupling constants defined on the bonds and re- 
spectively. The "staggered" arrangements of the bonds on a 
square lattice with periodic boundary conditions can be seen 
in Fig. [n^). The geometry of the ladder model results by a 
simple shift of every second dimer We define a = J'/ J as 
the parameter driving the phase transition. 

Simulations are performed with the directed loop variant 
1 19] of the stochastic series expansion (SSE) algorithm 1I20I1 
for lattice sizes L = 8 up to L 72 (in single cases L = 96) 
and inverse temperature up to /3 = 256. We checked that all 
quantities took on their ground-state values at the temperature 
simulated and we scaled /3 ~ L. Additional parallel temper- 
ing (PT) updates as well as multihistogram reweighting were 
performed to further optimize sample statistics and data anal- 
ysis. 

To probe the nature of the quantum phase transition, we cal- 
culate several well-known observables starting from the stag- 
gered magnetization (the Neel order parameter) with 



1 



N 



(2) 



and its Binder parameters Qi — {(m^)'^) / {\in^\)'^ and Q2 = 
{{mlY) / {{ml)'^)'^ . These quantities are complemented by 
the second-moment correlation length obtained from structure 
factors S as 



S'(7r,7r) 



S{-K,1T + 2'K/Ly) 



- 1. 



(3) 



with the obvious relation for the imaginary time correlation 
length ^T-. Lastly we determine the spin stiffness obtained 
from 



(4) 



where vo^ is the square of the difference of operator numbers 
S" and 5^5+ in x-direction. At a critical point the quan- 
tities Qi, Q2, S,y/L as well as psL are expected to cross for 
different lattice sizes L (under the assumption that the imagi- 
nary time exponent z = 1, in case of the spin stiffness). 



We first present QMC results for Q2 and iy/L in Fig. |2] 
where the crossing behavior becomes evident. A second-order 
phase transition is therefore very likely to happen. However, 
corrections to scaling terms are clearly present as the cross- 
ing points are not sharp but rather spread out for smaller lat- 
tice sizes. We exploit this fact by studying the scaling of the 
crossing points at lattice sizes L and 2L for the various quan- 
tities. In this way a bracketing of the critical coupling is 
obtained in Fig. |2jc) and we can easily read off a preliminary 
estimate as G [2.5190,2.5202]. This value is made more 
precise by fitting to a function a^{L, 2L) = + aL^^Z''"" 
yielding ac — 2.5198(3). All observables agree in this pic- 
ture, indicating a single phase transition and our estimate is in 
accordance to earUer quotes in the Uterature 1110051. 

Finite-size scaling: Having gained a fairly good estimate 
of the critical coupling ac we now turn to determining the 
critical exponent v. This is done by using the scaling ansatz 
for a second-order phase transition OL{t) = L^/^ go{tL^/'^), 
where A is the scaling exponent associated with the quantity 
O, t = ajoLc — 1 the reduced critical coupling and go a scal- 
ing function. For the quantities Q2 and ^yjL, it is clear that 
A = 0, which is indeed verified from the data in Fig. [J] In 
this work we follow Ref. H and take evident corrections to 
scaling explicitly into account by performing the data analy- 
sis according to a more general scaling ansatz 

OL{t) = L^'^il + cL-^)goitL^^'' + dL-'^'") , (5) 

where lo is the usual confluent correction exponent and a 
shift correction contribution. This way we can directly com- 
pare with a very detailed study recently performed on two bi- 
layer models favoring dimer formation which gave strong sup- 
port for 0(3) universality We perform our data analysis 
using this scaling ansatz in two ways. First, a Taylor expan- 
sion of go {x) up to fourth order in x is used in conjunction 
with multidimensional fitting. Second, we check this proce- 
dure by using a collapsing tool 12111 which makes direct use 
of multi-histogram reweighting. Both methods give consis- 
tent results for the critical coupling ratio as ac = 2.5196(2) 
and the critical exponent ly = 0.689(5) which is more than 4a 
smaller than the standard 0(3) value of j/ = 0.7112(5). 
We arrive at this result conclusively for all observables of this 
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Figure 3: (color online). Best collapse for the spin stiffness (upper 
data set) and tlie correlation length (lower data set) obtained from 
fitting data to the scaling ansatz l|5j. 



study. The error bar reflects checking for different window 
sizes, as well as trying different correction terms. In fact, we 
find that in most cases the oj correction is sufficient, i.e., in- 
clusion of (p terms does not change the estimate for i^. Fig- 
ure [3] contains a data collapse for all lattice sizes for quanti- 
ties and psL. The exponent z can best be estimated from 
the correlation length in imaginary time. As the result 
i^r = 0.687(5) is almost equal to the previous value we can 
conclude z ~ 1.01(1). 

Cross-checks and critical scaling: Since the discrepancy 
between ly and the standard 0(3) value is rather small we have 
performed multiple checks of the algorithm and our numer- 
ical procedure in different categories. First, we repeat the 
study for a different aspect ratio 2Ly x Ly {L^ = 2Ly) (and 
larger where the coiTelation lengths in the x and in the 
y direction are approximately equal, giving a consistent re- 
sult of I' = 0.688(5). Second, we carefully ran simulations 
on various dimerized models known to be described by expo- 
nents in the 3D Heisenberg universality class. These include 
the CaVO, the bilayer, and the ladder model discussed before. 
In all cases we arrive easily at Heisenberg universality. This 
is also true for the plaquette model of Fig. Hlb), being a model 
not previously investigated to high precision. There, we ob- 
tain a critical point of ac = 1.8228(4) and a critical exponent 
of f = 0.709(8) using exactly the same procedure (even at 
smaller lattice sizes of up to L = 48) 02611 . Those checks 
on known and hitherto less studied models indicate that the 
critical exponent ly is indeed smaller for the J- J' model. 

To further investigate the "mismatch" of the universality 
class we proceed with determining other critical exponents by 
studying the scaling at the quantum critical point ac- In this 
case the staggered magnetization scales as (|mg|) L"^/'^ 
and we can obtain the exponent 1] (as well as z) from {rril) 
and the staggered susceptibility from 



d — Z — 7] 



Xs 



(6) 



Using this approach, we obtain in Fig. |4ja) the estimates 
P/i/ = 0.515(4) for the bilayer, the ladder and the pla- 
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Figure 4: (color online). Scaling of (a) the staggered magnetization 
and the (b) the slope of ^y/L at the critical point. We compare the 
scaling of the J- J' model (jjp) against the plaquette (pl), the bilayer 
and the ladder models at the best known critical couplings. Both 
quantities indicate different critical exponents for the J- J' model. 



quette model at the known critical points (see Table and 
(ijv = 0.545(4) for the J-J' model, which should be com- 
pared to the 0(3) value of jijv = 0.518(1) The error 
bars on the data reflect uncertainties in ac and the straight line 
fits are all excellent and results are independent of different fit- 
ting windows. Our results are quoted for the five largest lattice 
sizes. To make the discrepancy in more apparent we have 
rescaled the original data to start at a common point in the plot. 
Second, we compute v again from the slope sq = AQ2 /da 
and = (l/L)d^y/da at the critical point which should 
scale with lattice size as Li/"^. Fi gure |4jb) shows this for 



Table I: Heisenberg dimer and quadrumer systems used for compar- 
ison of scaling at the critical point. 



Model 


Type 


Qfc Reference 


Bilayer 


Symm., dimer 


2.5220(1) [&] 


Kondo 


Symm., dimer 


1.3888(1) [&] 


CaVO 


Symm., plaquette 


0.939(2) [5] 


Plaquette 


Symm., plaquette 


1.8228(4) [2&] 


Ladder 


Unsymm., dimer 


1.909(1) [7,2^ 


J-J' 


Unsymm., dimer 


2.5196(2) this work, [2^ 
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Figure 5: (color online). The critical Binder parameter in dependence 
on the lattice size L for the (a) plaquette model, (b) ladder model, 
(c) bilayer, (d) Kondo lattice, (e) classical 0(3) model, and (f) J-J' 
model. Data for the CaVO lattice is not shown as they overlap with 
curve (c). 



in comparison for the J-J' , the ladder, and the plaquette 
model. Fits for the ladder and the plaquette model yield a 
common v ~ 0.709(6) while v = 0.689(5) is obtained for the 
J-J' model, in accordance with the previous analysis. Sim- 
ilarly the remaining exponents 77 and z are determined to be 
d-z~ri = 0.908(5) (J-J') as well as d - z - 77 = 0.971(2) 
(other models) with the obvious contrast. It is easily checked, 
that the scaling law 2(3 = (d + z — 2 + rj)i^ is satisfied within 
error bars for all cases. The exponent ?/ for the J-J' model is 
thus given by 77 = 0.09(1), which is considerably larger than 
the standard 0(3) value. 

Our findings are finally reinfored by comparing the Binder 
parameter at the best known critical points for the models of 
TableU It is evident from Fig.[5]that all cases apart from the J- 
J' model are in accordance with 0(3) behavior. To make this 
comparison even stronger we also include in Fig. |5]the value 
from Wolff cluster simulations of the ordinary 3D classical 
Heisenberg model ll29ll . 

Conclusion: In this Letter, we give comprehensive nu- 
merical evidence for an unconventional universality class of 
the J-J' model based on data collapsing analysis, scaling at 
criticality and by a comparison of the Binder parameter for six 
different dimerized models. This shows that there are nontriv- 
ial contributions to the quantum critical point changing the 
critical exponents. Those contributions are triggered by the 
special staggered arrangement of couplings. Our result chal- 
lenges the current understanding of quantum phase transitions 
in dimerized quantum spin systems and it will be interesting 
to see which exact theoretical mechanism accounts for the ob- 
served discrepancy. 
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